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tt-ura-.ru:. 


1 .  Overview  of  CMS-Wave 


■  Steady-state  (time-independent),  half-plane,  two- 
dimensional  spectral  transformation  solved  by 
finite-difference,  forward-marching  implicit 
scheme 

■  PC-based  efficient  model,  stand-alone  or  coupled 
to  CMS-Flow,  a  circulation  and  sediment  transport 
model,  through  the  SMS  interface 

■  Emphasis  on  wave-structure-land  interactions  for 
practical  coastal  engineering  projects 


Coastal  Modeling  System  Basics  Webinar 


3 


2.  Capabilities 


■  Wave  diffraction,  reflection  (forward  & 
backward),  breaking,  bottom  friction  dissipation 

■  Wind  input,  wave-current  interaction 

■  Wave  transmission  at  structures 

■  Wave  run-up,  overtopping,  overland  flow 

■  Variable  grids  with  nesting 

■  Nonlinear  wave-wave  interaction  &  infra-gravity 
waves 

■  “Fast  mode”  for  quick  calculations  &  prelim  runs 
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Structures 


CMS-Wave  and  STWAVE 


CMS-Wave  and  STWAVE  (half-plane)  Comparison 

Capability 

CMS-Wave 

STWAVE 

Spectrum  transformation 

Directional 

Directional 

Refraction  &  shoaling 

Represented 

Represented 

Depth-limited  wave  breaking 

Choice  among  four  formulas 

One  formula 

Roller 

Represented 

None 

Diffraction 

Theory 

Smoothing 

Reflection 

Represented 

None 

Transmission 

Formulas 

None 

Run-up  and  setup 

Theory 

None 

Wave-current  interaction 

Theory 

Theory 

Wave-wave  interaction 

Theory 

Semi-empirical 

Wind  input 

Theory 

Semi-empirical 

White  capping 

Theory 

Semi-empirical 

Bottom  friction 

Theory 

Theory 
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CMS-Wave  SMS  11.0  Interface 


S  11.0  64-bit  -  [untitled. sms] 
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3.  Governing  Equation 


Wave-Action  Balance  Equation  with  Diffraction 


8[(Cgx+u)A]  d[(c„+v)A]  d[ce9A]  k 


dx 


+ 


gy 


dy 


+ 


_go_ 

SO 


Diffraction  intensity  factor 

{(cccos20A)y -]-cc ■  cos20A\  +  Sin  +  S 


2a"  s  yy  2  8  w 


ap 


where  A  =  E/a  ,  wave-action  spectrum 
and  E  =  E(a,Q )  ,  wave  directional  spectrum. 


Note:  x  is  normal  to  the 
offshore  boundary; 
y  is  parallel  to  the 
offshore  boundary 


y 

A 


■-vv-v 

*/vv. 

•t-V-V 

■'.V.V. 

■■■■■■ 

vv 


x 


shore 


.y.J 

v-v:-:-: 

V;/./ 

'.i  ■ 

vvy 


Coastal  Modeling  System  Basics  Webinar 


7 


4.  Incident  Wave  Spectrum 


■  NDBC/NOAA  Ocean  Buoys 


■  CDIP  Coastal  Buoys 

■  Project  specific  measurements  (ADCP) 


■  Theoretical  spectra  (SMS) 
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SCO  (m2/Hz) 


Theoretical  Spectrum 


A  single  input  spectrum  applied  along  the  seaward 
boundary, 

e.g.,  a  JONSWAP  type: 


E  = 


ag 

<75 


<J, 


exp(-0.74  -^-)yaD((j,  6) 


u 


where 


D(  ff)  = 


2s  r(s/2+i) 
n  r(s + 1) 


coss(0-0o) 


for  1 6  -  0o  \<  k  /  2 


and  s  is  the  directional  spreading  parameter. 
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SCO  (m  2 /Hz) 


Idealized  Directional  Distribution 


-100  -80  -60  -40  -20  0  20  40  60  80  100 
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SMS10.1  Wave  Spectrum  Display 


63 
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Spectral  Viewer 


View:  j Polar  Energy: 

Spectral  Energy 
-  10.0 
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5.  Wave-Current  Interaction 


■  Solving  for  wave  number  k  in  dispersion  equation  with 
a  current: 


<7  =  y]gk tanh kh  +  ku  cos#+£vsin# 

■  Computing  wave  radiation  stresses 


sxx  =E[n(cos2e  +  \)-^\, 
Syy=E[n(  sin20  +  l)-|], 

A* 


~n  .  „„  ,  1  kh 

2? —sin 20,  where  n  =  — l - 

2  2  sinh  2  kh 
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6.  Jetty  Breakwater  Wave  Diffraction 

and  Reflection 


water  surface 


(b)  Elevation 


Coastal  Modeling  System  Basics  Webinar 


13 


Infra-gravity  Waves  at  Humboldt  Bay,  CA 
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Wave  Period  (sec)  Wave  Height  (m) 


tt-ura-.ru:. 


7.  Wind-Wave  Generation 
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Hs  (cm) 


Wave  Breaking  Formulas 
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Wave  Generation  in  Matagorda  Bay,  TX 
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Variable  Rectangular-Cell  Grids 


Variable-rectangular  cells 
Total  223  x  172  cells 


Square  (20  m  x  20  m)  cells 
Total  316  x  426  cells 
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CMS-Wave  on  Variable  Grids 


Variable-rectangular  cells 
Total  223  x  172  cells 


Square  (20  m  x  20  m)  cells 
Total  316  x  426  cells 
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Grid  Nesting 
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Regional  Wave  Generation 

Incident  Waves:  12.9  m,  13.8  sec,  from  S 


Max  Surge:  3.5  m  (Return  Period  =  50  yrs) 


Without  wind 


With  wind  (27  m/sec,  from  S) 
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8.  Wave  Run-up 


Wave  run-up:  rush  of  waves  up  a  slope  or 
structure 

Two-percent  run-up,  R2  :  the  vertical  up-rush 
level  exceeded  by  2-percent  of  the  larger 
run-up  height 

Ahrens  &  Titus  (1981),  Mase  &  Iwagaki  (1984) 
-  400  laboratory  experiments 
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Wave  Run-up  Calculation 


Total  run-up  R2  =  wave  setup  +  2%  exceedance  of  swash  level 


Wave  setup: 


o]l  = _ 1  ,dSa  |  dS  drj_  = _ L(^  +  ^) 

dx  pgh  8x  dy  '  dy  pgh  dx  dy 


Max  setup  (Guza  and  Thornton,  1981):  ^inax  =  0.17//0 

Total  runup  R2  (2%  exceedance)  =  2  7max  (Komar,  1998) 
Max  water  level  =  max  of  (  fj  +  Hs/2  ,  R2) 

*  Wave  setup  and  max  water  level  field  are  saved  in  setup.wav 
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Specify  Feature  Cells  in  SMS11.0 


11.0  64-bit  -  [untitled. sms] 


g  File  Edit  Display  Data  CM5-WAVE  Web 


-  I  n  x 
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Floating  Breakwater 


An  analytical  formula  of  the  transmission  coefficient 
for  a  rectangle  floating  breakwater  of  width  B  and 
Draft  D  (Macagno  1953): 


1+ 

kh  ^ 

kB  sinh 

2  n 

2 

2  cosh  k(h-D) 

J 
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Bottom-Mound  Breakwater 


Vertical  wall  breakwater  (Kondo  and  Sato,  1985): 

K  =  0.3  (1.5 — — ),  for  0<-^<1.25 

h;  hs 

Composite  or  rubble-mound  breakwater: 

K  =0.3  (1.1— -M,  for  0<— ^<0.75 

H:  H . 

s  s 

where  hc  is  the  crest  height  (above  mean  water  level) 
and  H s  is  the  incident  wave  height. 
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Idealized  Island  Example 


20  feature  cells 

input  depth 
=  10  m 

incident  wave: 

2  m,  6  sec, 

30  deg  oblique 
(gamma  =  4) 
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Idealized  Floating  Breakwater 


20  feature  cells 

Input  depth 
=  10  m 


incident  wave: 
2  m,  6  sec, 

30  deg  oblique 
(gamma  =  4) 

draft  =  2  m 
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Idealized  Platform 


® 


20  feature  cells 

input  depth 
=  10  m 


incident  wave: 
2  m,  6  sec, 

30  deg  oblique 
(gamma  =  4) 

platform  elev. 

=  1  m  (mwl) 
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Submerged  Platform 


20 

9  10  4-2 

10  10  4-2 

11  104-2 

12  10  4-2 

13  10  4-2 
9  11  4-2 

10  11  4-2 

11  11  4-2 

12  11  4-2 

13  11  4-2 
9  15  4-2 

10  15  4-2 

11  15  4-2 

12  15  4-2 

13  15  4-2 
9  16  4-2 

10  16  4-2 

11  164-2 

12  16  4-2 

13  16  4-2 


struct 


20  feature  cells 

input  depth 
=  10  m 


incident  wave: 
2  m,  6  sec, 

30  deg  oblique 
(gamma  =  4) 

platform  elev. 

=  -2  m  (mwl) 


Coastal  Modeling  System  Basics  Webinar 


30 


TrtJiMitsion  Coefficient ,  Kr^Hrf  H> 


Wave  Transmission  Experiment  (Goda,  2000) 


Transmission  coefficients  kt 

H j  =  1  m,  Tp  =  6  sec  (monochromatic  wave) 

/)  =  10  m,  d=5m,  B  =  80  m 

K  (m) 

CMS-Wave 

Equations 

Vertical 
wall  ^ 

Rubble 

mound 

Vertical 

wall 

Rubble 

mound 

-2.0 

1.02 

1.02 

-1.5 

1.03 

1.03 

-1.0 

0.78 

0.78 

-0.5 

0.63 

0.63 

0.0 

0.46 

0.34 

0.45 

0.33 

0.5 

0.27 

0.18 

0.30 

0.18 

1.0 

0.15 

0.04 

0.15 

0.03 

1.5 

0.10 

0.024 

2.0 

0.07 

0.018 

Wave  overtopping:  Surge  level  =  0.81  m  (3  ft) 
Hs  =  0.88  m,  Tp  =  10.1  sec  (Hughes,  2008) 


ERDC/CHL  TR-08-1 0 
by  Hughes (2008) 
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Water  surface,  m 
7.0 

6.5 

^  6.0 

5.5 
5.0 

4.5 
4.0 

3.5 
3.0 

—  2.5 
2.0 

1.5 
1.0 


Calculated  Wave  Overtopping  R127 
Surge  level  =1 .3  m,  Hs  =2.3  m,  Tp  =14  sec 


Coupled  CMS-Flow 
and  CMS-Wave 
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Calculated  Wave  Overtopping  Rate 


Case 

number 

Surge 
level  (m) 

Wave 
height  (m) 

Wave  peak 
period  (sec) 

Overtopping  rate  (m2/sec) 

Measured 

CMS-Flow 

CMS-Wave 

R128 

0.29 

0.27 

0.28* 

0.29 

0.82 

6.1 

0.38 

0.38 

0.39 

R109 

0.29 

0.26 

0.28* 

0.29 

2.48 

13.7 

0.70 

0.85 

0.92 

R121 

1.3 

2.55 

2.57* 

1.3 

2.30 

6.1 

2.67 

2.93 

2.76 

R127 

1.3 

2.54 

2.57* 

1.3 

2.31 

14.4 

2.84 

2.98 

2.81 

*  Calibration  □  With  wave  overtopping 
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Muddy  Bottom 


Wave  dissipation  by  damping  (Lamb,  1932): 

sdP=-^vk+v^klE 


where  V, 


is  the  kinematic  viscosity  of  sea  water, 


and  vt  is  the  turbulent  eddy  viscosity: 


=v 


Hl 

t, breaking  fa 
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Wave  Height  Hs,  m  Wave  Height  Hs, 


Louisiana  Muddy  Coast  Simulation 


Dsptft.  m  (MSL) 

—  1DDD 

—  500 

—  400 

—  300 
-  200 

—  too 

—  0 


J 

Primarily  sand 


42041 

— - 


Mud  area 


January  2004 
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® 

— 


CMS-Wave  Fast  Mode 

(Simplified  Formulation) 

~ 


•  Fast  mode  uses  5  to  7  directional  bins  with  spectral  calculations 
(Standard  runs  with  35  directional  bins) 

•  Ideal  for  quick  applications,  prelim  runs,  time-pressing  project 


Standard  run 


Fast  mode 
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Nonlinear  Wave-Wave  Interaction 


Governing  Equation: 


DA 

Dt 


^diffraction  ^ in  ^ dp  , nl 


where  Snl  is  the  nonlinear 
wave-wave  interaction  term 


dB  B 

Anisotropic  Snl:  Snl=a(a)— — I- Z>(cr) -^—5-  (Jenkins &  Phillips,  2001) 


da 


se- 


1  a 

where  a- — —  [1 H-  (2/2  —  l)2  cosh2M]-l,  b-  — 


2  n 


no 


3 _ 5 


and  B  =  k o 


n 


(2  TcYg 


[&)4Ef 


2  LVJ 
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(m  sec) 


® 


Exact  and  Calculated  Snl(f  ) 


x  10' 


o 

ID 

</) 


(O' 


0.8 

0.6 

0.4 

0.2 

0 

-0.2 

-0.4 

-0.6 

-0.8 

-1 


Frequency  (Hz) 


r=5 


Frequency  (Hz) 
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Direction  (deg)  Direction  (deg) 
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Spectral  Evolution  and 
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Nonlinear  Wave  Effect 
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9.  Coupling  with  CMS-Flow 


Breaching  at  Jetty,  Simulation  at  Matagorda  Ship  Channel,  TX 


Morphological 

Change 


\ 


E  South 
Jetty 

Breaching 
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MSC  Jetty  Wave  Run-up  &  Breaching 
Cat  3  Hurricane  (50-Yr  Life-Cycle) 


N.  Jetty 
breaching 


S.  Jetty 
breaching 


After  12-hr  simulation 


Initial  bathymetry 


•  Peak  storm  surge  level  reaches  3.5  m  between  Hrs  4  and  8 

•  Incident  offshore  wave  is  7.6  m,  14.3  sec,  from  south 
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MSC  Jetty  Wave  Run-up  &  Breaching 
Cat  3  Hurricane  (50-Yr  Life-Cycle) 


Storm  surge  over  the  initial  bathymetry  South  Jetty  breach  in  12-hr  simulation 


•  Peak  storm  surge  level  reaches  3.5  m  between  Hrs  4  and  8 

•  Incident  offshore  wave  is  7.6  m,  14.3  sec,  from  south 
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Calculated  30-day  Morphology  Change 
Tombolo  Development 


Depth,  m 


CMS 

Steering  Interval 
=  4  hr 


Grain  Size 
=  0.18  mm 

Hydro  time  step 
=  0.25  sec 

Transport  and 
morphology 
calc  time  step 
=  9  sec 


□  DD:QQ:00 
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10.  Future  Development 


■  Telescoping  grids 

■  Dynamic  memory 

■  Full-plane  transformation 
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Conclusions 


■  CMS-Wave  designed  for  wave-structure-land 
interactions  for  inlet  and  nearshore  applications 

■  Coastal  inlet-specific  processes  represented 

■  Emphasis  on  computational  speed  and  SMS 
integration  for  PC  users 

■  Coupled  to  CMS-Flow  for  sediment  transport 
and  morphology  change 
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